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ABSTRACT 

Using recent results of Avdeev et al. and an expansion for fj, t /rht(Mt) (M t is the 
pole mass and /i t = rh t (fi t ) ), it is shown that when Ap is expressed in terms of 
rh t 2 (M t ), the QCD correction is only (2 — 3) x 10~ 3 in the NLO approximation. As 
a consequence, in terms of M 2 the correction to Ap is almost entirely contained in 
rh t 2 (M t ) I , a pure QCD effect. The latter is studied using various optimization 
procedures, and the results compared with the expansion proposed by Avdeev 
et al.. Implications for electroweak physics are discussed. Threshold effects are 
analyzed on the basis of a simple sum rule. 



The question of the QCD corrections to Ap has been the subject of several recent 
studies [1-4]. In particular, Avdeev, Fleischer, Mikhailov, and Tarasov carried out a 
complete three-loop calculation of 0(aa 2 ). Their result, obtained in the limit m b — > 0, 
can be expressed to good accuracy in the form 



(A,)/ = ^[l + « 



QCD\ 
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where / denotes the fermionic contributions, M t is the pole mass, and 
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Sqcd = -2.860^^ - 10.55 J (2) 

is the QCD correction. As an illustration, for M t = 200GeV, Eq.(fJ) gives Sqcd = 
-0.0961 - 0.0119 = -0.1080. Thus, the 0(a 2 s ) term in Eq.(@) is reasonably small 
and leads to an enhancement of 0.0119/0.0961 = 12.4% of the leading QCD result. 
Moreover, if a s (M t ) in Eq.(^) is expressed in terms of a s (fi) and the resulting series is 
truncated in 0(a 2 (p)), the //-dependence of Sqcd{h) is relatively mild. For example, 
for M t = 200GeV and 0.1 < /i/M t < 1, we find a variation < 5.8 x 10~ 3 , which 
amounts to 5.4% of the total QCD correction or 49% of the 0(a 2 ) term. 

In spite of these facts, it is worth noting that the expansion in Eq.(||]) involves 
large and increasing coefficients, a feature that is frequently an indication of significant 
higher order terms. Furthermore, the arguments of Ref. [[J suggest that there are at 
least two scales in this problem: one, of 0(M t ), associated with the correction to the 
electroweak amplitude, and another one, much smaller, related to contributions to the 
pole mass M t involving small gluon momenta. It should also be observed that the /i- 
sensitivity of the truncated series depends on the chosen interval, and becomes sharper 
for fi/M t < 0.1. For these reasons, it is a good idea to find alternative expressions for 
Sqcd that separate the two scales, and at the same time involve terms of 0(a 2 ) with 
small coefficients. A simple way of implementing this idea has been outlined in Ref. |J. 
One expresses first (Ap) f in terms of rht(M t ), the running mass evaluated at the pole 
mass, and then relates m t {M t ) to M t by optimizing the expansion of M t /rht(M t ), 
which is known through 0(a 2 s ) 0. As shown below, the arguments of Ref. p| can be 
significantly refined using the new results of Ref. Q . 

Calling p t the solution of m t {[j) = fi and using Eq.(19) of Ref.||, (Ap)j can be 
written in the form 

6%B D = -0.19325 ^M + o.0711l(^^] . (4) 
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We see that the convergence pattern of Eq.(4) is very nice, with very small and decreas- 
ing coefficients and alternating signs. For this reason we will assume that the terms of 
0(«g) and higher are negligible and evaluate Eq.(|) with \i t — > M t as this introduces 
only a small change of 0{ofy. This argument is further strengthened in Appendix 
A, where we discuss the potential contribution of threshold effects to perturbative 
expansions, such as Eq.(|4]), on the basis of a simple sum rule. 

In order to express (Ap) / in terms of m t (M t ), we use the following next-to-leading 
order (NLO) expansion, derived in Appendix B: 

* l+ U^)\ [36M _ 2A5ni] (^m\\ (5) 



m t (M t ) 3 V 7T J \ TT 

where n/ = 5 is the number of light flavors. Writing 



3G^m t 2 (M t ) 

(Ap)/ = — |H=^ — [1 + Aqcd], (6) 



it is possible to combine Eq.(3-5) in a single expansion 
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Aqcd = -0.19325 ^ + C { ] . (7) 

7r \ n J 

In fact, Eq.(0) was proposed in Ref.0 before the results of Ref. Q were known, and 
it was argued, on the basis of convergence assumptions, that \C\ < 6. From Eqs.(3- 
5) we see that C = 5.40, consistent with the arguments of Ref. 0. However, with 
C = 5.40 the two terms in Eq.(|7|) nearly cancel and its convergence properties become 
problematic. In order to evaluate Aq C d, h is better to employ the separate expansions 
in Eqs.(3,4) and Eq.(|5|). In the latter we retain the relatively large 0((a s /7r) 3 ) term in 
order to control the scale of the leading contribution. Applying the Brodsky-Lepage- 
Mackenzie (BLM) optimization || to Eq.(^[), we have 

fit 1 | 8 f a s (0.252M t ) \ 2 ^ ^ / q s (0.252M t 



rht(M t ) 3 V vr 
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For M t = 200GeV, Eqs.(|) and (|) give p t /m t (Mt) = 1.00392 and 1.00425, respectively 
(in this paper a s (p) is evaluated with a 3- loop /^-function for rif = 5 light flavors, 
normalized such that a s (mz) = 0.118, with rriz = 91.19GeV). We have verified that 
the evaluations of Ht/rht(M t ) by the priciple of minimal sensitivity (PMS)0 and the 
method of fastest apparent convergence (FAC) || are very close to the BLM result (the 
difference is < 5 x 10~ 6 ). The correction Aqcd i n Eq.(||) is obtained from the relation 



1 + Aqcd - ( —jjj- 
\m t {M t 



1 +tf MS 
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where Sqq D and /i t /m t (M 4 ) are evaluated via Eq.(£|) (with /i t — > M t ) and Eq.(||), 
respectively. 

Table 1 shows that Sqq D and Aqcd are indeed small corrections for M t > 130GeV. 
In particular, Aq C d = (2 — 3) x 10~ 3 , depending on M t . This conforms with the idea 
that m t (M t ) is a good expansion parameter in the sense that the associated QCD 
corrections are small 0. Furthermore, in the (Ap)j case this occurs to a rather re- 
markable degree. However, it is not a good approximation to retain only the 0(a s ) 
correction in Eq.(6), as it is nearly cancelled by higher order contributions. On the 
other hand, the expansion in Eq.@ is very well behaved. In this sense, fit is a bet- 
ter expansion parameter than rh t {M t ), which may reflect the fact that it makes no 
reference to the pole mass M t .[] 

Following the strategy outlined in Ref.||, Sqcd in Eq.(l) can be obtained from the 
relation 




1 + Sqcd = [-^f 1 ) [1 + A QCd}- (10) 
For M t /rht{M t ) we have the well-known expansion || 



M t 4a s (M t ) , 1 . (a s {M i 
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1 + — il_ IL + (16.11 - 1.04w f ) v T/ , (11) 
m t {M t ) 3 7T 1 f '\ ti ) 1 K ' 

where the masses of the light quarks have been neglected. As M t /rh t (M t ) plays a 
central role in our final result, we record the expressions obtained by applying the 
three optimization methods to the r.h.s. of Eq . (|TT|) 

«. 1 + i^l_ lm (M 2 (BLMl (12) 



rht(Mt) 3 7T \ 71 
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1 + — i^^_o, 8 4 1 (PATS'), (13) 



rht(M t ) 3 7r y 7r 
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m t (M t ) 3 7T 

where //* = 0.0963Mt@, //** = 0.1004Mi, /i*** = 0.1183M t . We note that Eqs.(12- 
14) involve coefficients of 0(1) and similar scales. For M t = 200GeV, Eqs. (12-14) 
give M t /rh t (M t ) = 1.06305, 1.06304, and 1.06299, respectively. For M t = 174GeV, 
the corresponding values are 1.06479, 1.06478, 1.06472. Finally, for M t = 130GeV, we 

1 The usefulness of /i t as an expansion parameter has also been emphasized in Ref.Q. 
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have 1.06877, 1.06876, 1.06869. Thus, the three approaches give remarkably consistent 
results. In contrast, the expansion in Eq . (|TT|) , which involves a large second order 
coefficient, gives, for M t = (200, 174, 130)GeV, 1.0571, 1.0584, 1.0614, respectively, 
which are (0.6 — 0.7)% smaller. One must conclude that the coefficients of the unknown 
terms of 0((a s /ir) 3 ) and higher in Eq . (|TTD and or Eqs. (12-14) are large. For instance, if 
the truncated BLM expansion in Eq.(|T2"D were exact, the coefficient of the (a s (M t )/7r) 3 
and (a s (M t )/ V) 4 terms in Eq. ([TT|) would be ~ 104 and 1,041, respectively.^ In the 
following we employ the optimized expression for M t /rht(M t ), which, for defmiteness, 
we identify with Eq.(|i~2j). 

Table 2 displays the values of Sqcd obtained from Eq. ([[(]), using Eq.(|i~2"D and our 
previous determination of Aqcd (Eq.©), an d compares them with those derived from 
Eq.(|2|). In order to show the effect of the higher order contributions (H.O.C), we also 
exhibit the fractional enhancement of the total QCD correction over the conventional 
0(a s ) result (first term in Eq.(Q)). 

As pointed out by Kniehl pTT | , the various evaluations of QCD corrections can be 



conveniently interpreted as effective re-scalings of the 0(a s ) result. In fact, we find 
that the two evaluations of Sqcd in Table 2 can be very well represented by the simple 
formulas 

Sqod = . 2M ^M^1 lEq , m (15) 

7T 

a s (0.323M 4 ) . 
Sqcd = -2.86^ ^ Sg.(10). (16) 

7T 

The effective scales have been chosen so that Eqs. ( ]T5|JT6| ) reproduce the values in Table 
2 with an error of at most 1 or 2 times 10 -4 . 

From Table 2 we see that, for 130GeV< M t <220GeV, the evaluation of |<5qod| 
based on Eq.(|H]) is (5.1—6.5) x 10~ 3 larger than the results from Eq.(|2|). This difference 
amounts to m 5% of the total QCD correction. On the other hand, the last two columns 
in Table 2 show that the H.O.C. are ~ 45% larger in the evaluation based on Eq. (|T0"|) . 
As the two calculations coincide through terms of 0((a s (M t )/n) 2 ), one must conclude 
that the coefficients of the 0((a s /7c) 3 ) and higher terms in Eq.(fJ) and or Eq.flTOD are 
large. For instance, if Eq.fllOl) (with Aqcd and M t /m t (M t ) evaluated via Eqs.© and 
( [L2|)) were exact, the coefficient of the neglected (a s (M t )/ii) 3 term in Eq.(|]) would 
be ~ —93.4. Of course, it could happen that the exact answer lies between the two 
evaluations, in which case the 0(a 3 ) terms would be smaller. However, because of its 
pattern of large and increasing coefficients, we think that this problem is more likely to 
occur in Eq.(^). It is also useful to note that if one applies the optimization procedures 
to Eq.(|^), the results become closer to those from Eq. ([T0|) . Specifically, we have 

2 For a recent application of the PMS and FAC approaches to estimate higher order coefficients 
in other cases, see Ref. 110f| . 
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Sqcd = -2.86 
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(Eq.(2); BML), 



(17) 



^ = -2.86^^ + 1.80 



7T 



7T 



(Eg. (2); PMS), 



Sqcd = -2.86 



a s (0.382Mi) 



TT 



(£g.(2); FACT). 



(19) 



For M t = 200GeV, Eqs. give 0.1085, 0.11044, 0.11037, respectively. We see that 
the difference between Eqs.( |l8j,|T9j) and Eq. fllOD is ~ 45% smaller than that between 
Eq.(Q) and Eq. fllCf) . It has also been pointed out that in most NLO QCD calculations, 
it is a good approximation to retain only the leading term, evaluated at the BLM scale 
[2,7]. From Eq.fllTD we see that this is not the case for (Ap)/, as the coefficient of the 
residual 0((a s /7i) 2 ) term, 9.99, is quite large. 

There is a simple but important point that should be stressed: once it is recognized 
that 5qcd is almost entirely contained in the first factor of Eq.fllOD (and this follows 
from Eq.(4) and (8), which are quite robust), it becomes clear that its magnitude 
and precision are largely controlled by the value of (Mt/rht(M t )) 2 and the accuracy 
within which it can be calculated. We emphasize that this is a pure QCD, rather 
than electroweak effect. We have also pointed out that there is a significant difference 
between the evaluations of this correction via Eq . (|i~l~D or by means of Eqs.(]T2|-[l4|). In 
this paper we have taken the point of view that (M t /m t (M t )) 2 can be best evaluated at 
present by means of the optimized expansions of Eqs. fll2| - |r4"D . In particular, we have 
noted that the coefficients in these equations are 0(1) and that the values derived 
from the three optimization procedures remarkably consistent. In fact, it would be 
very interesting to check this approach by an explicit evaluation of the unknown third 
order term in Eq.([TT|). 

As shown in Table 2, the formulation of this paper, based in Eqs. (|],|ll],|l!|) and 
summarized in the effective formula of Eq. (|T6|) , leads, for M t = 200GeV, to an en- 
hancement of 18.0% relative to the conventional 0(a s ) calculation. This modification 
corresponds to an additional contribution —2.2 x 10~ 4 to (Ap)j and +0.76 x 10~ 3 
to Ar, and shifts the predicted values of M t and rri\y by +1.9GeV and — 13MeV, re- 
spectively. It can be compared with an enhancement of 12.4% from Eq.(Q), or 14.9% 
from the optimized expressions of Eq . fll8|) or Eq . (|T9|) . Although the H.O.C. obtained 
from Eq.(|) and Eq.([K| are significantly different, the effects of this difference on 
electroweak physics, namely variations of 0.6GeV in M t and 4MeV in m\y, are rather 
small. As mentioned before, it is possible that the exact answer lies between the two 
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calculations. Although we prefer the evaluation based on Eq.(lO), the difference with 
Eq. (2) is interesting because it may be used as an estimate of the theoretical error due 
to uncalculated higher order terms. We recall that there is at present an additional 
uncertainty of ~ 5% due to the ±0.006 error in a s (m^)0. Following the discussion 
of Ref. 0, it is also easy to see how the approach in this paper can be applied to other 
electroweak amplitudes proportional to M 2 , such as those present in Z° — > bb. 

Finally, it is instructive to compare the above results with some of the estimates 
made before the complete 0(a 2 ) contribution became known. Although conceptually 
very different, the dispersive approach as implemented in the last paper of Ref. Jy], 
gives, for 130GeV< M t <220GeV, results very close numerically to those of Ref. [0]. 
The latter approximates the relevant Feynman diagrams by the 0(a s ) contribution 
evaluated at the BLM scale (first term of Eq.(17)). For M t = 200GeV, both estimates 
led approximately to a 34% enhancement of the 0(a s ) contribution, a result of the 
same sign but clearly larger than the conclusions from either Eq.(2) or Eq.(10). On the 
other hand, out previous estimate [(26 ± 6)%] ||, is roughly consistent with Eq.(10), 
but larger than the values from Eq.(2). 
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APPENDIX A 

We briefly discuss threshold effects on the basis of a simple sum rule that follows as a 
corollary from the arguments of Ref. fl2| . In that work the operator product expansion 
[ p"3|1 is applied to investigate the validity of the dispersion relations [D.R.] for vacuum 
polarization functions previously proposed by Chang et al. and Kniehl and Sirlin [p]], 
to all orders in perturbation theory. We call Yl^' l f(q,m 1 ,m 2 ) the vacuum polarization 
tensors associated with the vector and axial vector currents attached to pairs of quarks, 
such as ti, tb, bb, endowed with masses mi and m 2 . Decomposing 

n ^ A (<?, m u m 2 ) = Tl v ' A (q 2 , nix, m 2 )g iiv + \ v ' A (q 2 , m h m 2 )q il q vi (A.l) 

we have the sum rule 

o {^ V ( S '> m *> m i) + ^ A ( S '> m i)) 




oo 

ds' Im 
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A y (V, m 1 ,m 2 ) - A A (s', m 1 ,m 2 ) 



0. 



(A.2) 



In order to derive Eq.(A.2), one invokes analiticity and applies Cauchy's theorem to the 
appropriate combination of A's over a contour with straight lines just above and below 
the positive real axis, closed by a large circle. The arguments of Ref. ]I2] and the results 
of Ref. then show that the contribution of the large circle vanishes in the limit of in- 
finite radius. Identifying rui = m t , m2 = mj, there are large and equal non-relativistic 
contributions from the threshold region to ImX v (s' ,m t ,m t ) and Im\ A (s' ,m t ,m t ) 
(threshold contributions to the other A's are suppressed by reduced mass effects). 
However, the appropriate D.R. |IJ involve the integrals J°° ds'Im\ A (s' ,m t ,m t ) or 
J°° ds'Im\ v \s' ,m t ,m t ) over all values of s'. It is precisely from the contributions of 
such integrals in the D.R. that large threshold effects may potentially arise. On the 
other hand, Eq.(A.2) tells us that J°° ds'Im[\ v (s 1 , m t , m t ) + X A (s', m t , m t )] can be ex- 
pressed as a linear combination of integrals of the spectral functions ImX^'(s', m t , m b ) 
and Im\V' (s 1 , m&, m&) (j = V,A) which do not have significant threshold contribu- 
tions and involve different channels. The most sensible interpretation of this result is 
that the large threshold contributions to ImX v \s' ,m t ,m t ) and X A (s' ,m t ,m t ) cancel 
against other contributions when the integrals over all s' values are considered. As 
it is known that such integrals can be expanded perturbatively in powers of a s , this 
cancellation must occur order by order in pertubation theory. In particular, this ap- 
plies to the expansion in Eq.(4). Similar conclusions were reached in Ref.0 on the 
basis of different arguments. In summary, although the sum rule in Eq.(A.2) does not 
provide a rigorous proof of the cancellation of large threshold effects in electroweak 
amplitudes, it gives strong support to such conclusion when all contributions of a given 
order are included. On the other hand, it does not reveal the detailed mechanism of 
the cancellation, an interesting problem which hopefully can be clarified in the future. 



APPENDIX B 



In order to derive Eq.(5), we start with the well-known expression (TJj 



rutin) = m* -fa 



oisin) 
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1 + CL\ r 0-2 



+ 



71 



71 



(b.i; 



where ml is the RG invariant mass, (3\ = -j- — j, d = —2/ Pi, ai = iP 2 j Pi)\2j Pi — 
I2/P2}, 72 = — ^ (the other quantities need not be specified for our purposes). 
Setting fi = fit = m tint) m Eq.(B.l) and dividing by m t (M 4 ) we have 



rfitiMt 



l + ai^ + a 2 (^) 2 + 



l + ai ^(Md + a2 f^M t ) 



(B.2) 
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In the numerators we expand a s (/x) in powers of a s (M t ). The first factor gives 




1 + d^lJ^f 

71 \M t . 



Pi + P: 




(B-3) 



where we have used the fact that ln(// t /M t ) = 0(a s ). Similarly, we see that the second 
factor in Eq.(B.2) becomes 1 + ai/?i(a s (M t )/7r) 2 ln(/i 4 /M t ) + 0(aj). Combining the 
two factors, noting that (3\d = —2, fad + aifa = —72, and splitting the logarithm we 
have 



m t (M t 



1 - 



a 8 {M t . 



71 



In 



ft 



m t (M t 



+ In 



'rht(M t 



2 + 72 



a s {M t 
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(B.4) 



Inserting /j, t /m t (M t ) = 1+E«=i x n (a s (M t ) /i\) n in both sides of Eq.(B.4), using Eq.(ll) 
in the second logarithm and matching coefficients, we obtain x\ = 0, x<i = 8/3, 
x 3 = 2(16.11 - 1.04/2/ - 8/9) + 472/3 - 2^2 = 36.33 - 2.45n /; which leads to Eq.(5). 
This derivation assumes that rh t {^) evolves according to the nf light quarks. The 
difference in 72 and £3 when one employs 6 active quarks, as in Ref. is very small. 
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M t (GeV) 


io 3 fe 




130 


-6.80 


2.98 


150 


-6.67 


2.65 


174 


-6.53 


2.33 


200 


-6.41 


2.06 


220 


-6.33 


1.88 



Table 1. The corrections Sq- cd and Aqcd- The first one is given by Eq.(4) with /j, t — > 
M t , while the second is obtained from Eq.(9), with /i t /rh t (M t ) evaluated according to 
Eq.(8) (a s (m z ) = 0.118 is employed). 



M t 
(GeV) 


Sqcd 
(Eq.(2)) 


8qcd 
(Eq.(10)) 


H.O.C. 

(Eq.(2))% 


H.O.C. 
(Eq.(10))% 


130 


-0.1154 


-0.1219 


13.2 


19.5 


150 


-0.1128 


-0.1189 


12.9 


19.0 


174 


-0.1102 


-0.1159 


12.6 


18.4 


200 


-0.1080 


-0.1133 


12.4 


18.0 


220 


-0.1064 


-0.1115 


12.2 


17.6 



Table 2. Comparison of two determinations of Sqcd- The second column is based 
on Eq.(2) [4]. The third column is based on Eq.(10) with M t /rht(M t ) obtained from 
Eq.(12) and Aqcd evaluated according to Table 1. The fourth and fifth columns give 
the fractional enhancement over the conventional 0(a s ) result (first term of Eq.(2)) 
due to the inclusion of higher order contributions (H.O.C.). 
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